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We report the detection of a molecular outflow driven by the brown dwarf binary FU Tau. Using the IRAM 30 m telescope we 
observed the 12 CO(2-l) (CO) emission in the vicinity of FU Tau and detected a bipolar outflow by examining the wings of the CO(2- 
1) line as we moved away from the source position. An integrated map of the wing emission between 3 kms~' and 5 kms~' reveals a 
blue-shifted lobe at a position of ~ 20 " from the FU Tau system and at a position angle of ~ 20°. The beam size of the observations is 
11" hence it is not possible to distinguish between the two components of the FU Tau binary. However as optical forbidden emission, 
a strong tracer of the shocks caused by outflow activity, has been detected in the spectrum of FU Tau A we assume this component to 
be the driving source of the molecular outflow. 

We estimate the mass and mass outflow rate of the outflow at 4 x 10~ 6 M G and 6 x 10~ 10 M Q /yr respectively. These results agree well 
with previous estimates for BD molecular outflows. FU Tau A is now the third BD found to be associated with molecular outflow 
activity and this discovery adds to the already extensive list of the interesting properties of FU Tau. 
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^ 1. Introduction 

i Young brown dwarfs (BDs) occupy the mass regime between 
O stars and planets and are therefore significant to any theory de- 
+2 scribing activity in star forming regions. Thus they have become 
the subjects of increased scrutiny in recent years (Luhman 2012). 
i Their formation mechanism is at present much debated and in- 
deed it has been postulated that they may form by more than one 
t-H mechanism (Whitworth et al. 2006). The simplest idea is that 
they form in the same manner as low mass stars i.e. through the 
gravitational collapse of substellar mass cores (Padoan & Nord- 
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lund 2004). These cores occur directly by the process of turbu- 



^? lent fragmentation. In this scenario, BDs are just scaled-down 
. versions of low mass stars. Detailed studies of the circumstel- 
1— 1 lar environments of young BDs provide critical constraints to 
different formation mechanisms and are needed to identify the 
dominant mechanism. In particular, if BDs form like low mass 
. . stars we expect their accretion/outflow properties to be analo- 
gous. As a low mass star forms it displays a series of ubiquitous 
K^S observational properties, for example accretion disks, outflows, 
rN excess emission in the near-infrared and visual absorption. 

The observational evidences gathered to date in various 
wavelength domains indicate that young BDs show accretion 
and ejection behavior similar to low mass stars. For example 
they demonstrate T Tauri-like accretion (Jayawardhana et al. 
2003; Natta et al. 2004; Monin et al. 2010; Rigliaco et al. 2011) 
and both optical and molecular outflows, driven by BDs, have 
been detected. ISO-Oph 102 is a good example. It is an accretor 
with an observed accretion disk (Natta et al. 2002, 2004) and re- 
cent ALMA observations have detected millimetre sized grains 
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* Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



in its disk (Ricci et al. 2012). Its optical jet was discovered by 
Whelan et al. (2005) through the spectro-astrometric analysis of 
the [O]/l6300 emission line. Forbidden emission lines (FELs) 
like [OI]/16300 are important coolants in shocks and therefore 
good tracers of jets. Traditionally jets from classical T Tauri stars 
(CLTSs) are investigated by studying their FEL regions. Phan- 
Bao et al. (2008) also detected a CO molecular outflow driven by 
ISO-Oph 102. The orientation of the blue and red lobes agreed 
with the optical observations. 

The question of outflow activity in BDs is an important one, 
as a sufficiently efficient outflow activity could provide an ex- 
planation as to why the central object mass does not reach the H 
burning limit (Bacciotti et al. 201 1; Whelan et al. 2009; Machida 
et al. 2009). Molecular outflows are an important large-scale 
expression of jet launching. Indeed, molecular outflows were 
one of the first observational manifestations of this process to 
be studied (Reipurth & Bally 2001; Bachiller 1996). While gi- 
ant Herbig-Haro (HH) flows are optically visible and composed 
of many HH objects, each group representing different episodes 
of mass ejection, molecular outflows begin when the powerful 
bipolar jets accelerate and drive outwards the molecular gas in 
the vicinity of their parent star. Although it is accepted that 
they are powered by the primary jet from the protostar, the exact 
way in which the jet interacts with the molecular material is still 
uncertain (Cabrit et al. 1997; Downes & Ray 1999; Downes & 
Cabrit 2007). Molecular outflows are primarily detected in the 
CO molecule and thus millimeter observations have dominated 
the search for them. These outflows are mainly detected from 
Class and I low mass stars which are still embedded in their 
natal material. Observations of molecular outflows driven by the 
more evolved Class II CTTSs are much rarer (Cabrit et al. 201 1). 

As of today, only two detections of molecular outflows from 
optically visible young BDs have been made so far (Phan-Bao 
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et al. 2008, 201 1), although it is postulated that due to the colder 
environment of BDs, molecular outflows may be more common 
than in CTTS. We have conducted a survey of young BDs with 
the IRAM 30m telescope in the 12 CO(1-0) and 12 CO(2-l) to test 
this hypothesis (Whelan et al. 2013, in prep.). Their approach 
is to target BDs known to be accreting and which also show ev- 
idence of outflow activity primarily in form of FELs, in a mass 
range of 0.02 M Q to 0.13 M Q including a few very low mass stars 
(VLMSs). In this letter we report the detection of a remarkable 
molecular outflow in FU Tau, as part of our IRAM survey. 

FU Tau (04 ft 23 m 35 s 4, +25°03'03"05) is a BD-BD binary 
with a projected angular separation of 5".l or 800 AU at the 
distance to Taurus and a position angle (PA) of ~ 145° (Luh- 
man et al. 2009). Its membership of the Taurus molecular cloud 
has been known for some time (Jones & Herbig 1979) and it 
is situated in a relatively isolated region of the cloud. Luhman 
et al. (2009) give the spectral type of FU Tau A at M7.25, corre- 
sponding to a mass of 50 Mj up and the spectral type and mass of 
the companion at M9.5 and 15 M Jup respectively. The wide na- 
ture of the FU Tau binary challenges models which suggest that 
BDs form when their accretion is halted due to ejection from 
their natal clouds, since the system appears to have formed irre- 
spectively of dynamical interaction with nearby stars. A further 
intriguing property of FU Tau A is its over-luminosity with re- 
spect to other members of the Taurus star-forming region of the 
same spectral type (Luhman et al. 2009; Scholz et al. 201 1). The 
spectral energy distributions (SEDs) of both components show 
excess emission indicating the presence of circumstellar disks 
and their disks are classified as being Class II by Luhman et al. 
(2009). Furthermore, optical spectra clearly show that accretion 
is on-going in FU Tau A. Stelzer et al. (2010) estimated the mass 
accretion rate from both the Ha and He I (/15876) lines with Mhu 
= 3.5 x 10~ 10 M Q yr" 1 and M HeI = 7.5 x 10~ 10 M Q yr" 1 . Evi- 
dence of outflow activity comes from the detection of the [O 
l]A 5577 and [O l]A 6300 forbidden lines in an optical spectrum 
of FU Tau A (Stelzer et al. 2010). All of these facts combined 
show that the FU Tau system is probably a rarity amongst BDs 
and thus it is of considerable interest to test models describing 
the formation and evolution of BDs. For this reason we chose 
to publish the discovery of its molecular outflow separately from 
the global presentation of our IRAM survey. 

2. Observations and Data Reduction 

Observations of the CO(2-l) line emission were carried out at the 
IRAM 30m telescope on July 16th- 18th 2011 using the EMIR 
receivers at 1.3mm. 

In a first step, deep integrations were performed towards the 
protostar and at a reference position located 20" away. In a 
second step, in case of significant variations of the CO emis- 
sion between both positions, i.e. beyond the 3cr intensity level, 
more extended mapping at 12" sampling was performed. The 
CO emission map detected toward FU Tau is displayed in Fig 1 . 

Observations were carried out in Frequency Switch mode us- 
ing a throw of 14.3 MHz at 1.3mm, with a phase time of 0.2 sec- 
ond. An autocorrelator providing us with a spectral resolution of 
40 kHz was used as spectrometer. The weather conditions were 
rather good and stable, with system temperatures r sys varying 
between 200 and 400 K. Each position was observed so to reach 
a final rms of about 40 mK unless exception per velocity interval 
of O.lkm.s -1 , after averaging both polarizations. 

Pointing was checked every 1.5 to 2 hours and was found to 
be very stable, with pointing offsets corrections less than 3". The 
telescope parameters are adopted from the IRAM webpage. At 



the frequency of the CO(2-l) line, the main-beam efficiency of 
the telescope is 0.59 and the half-power beamwidth is 11". The 
intensities of the measurements are expressed in units of main- 
beam brightness temperature T m \,. 

The data were reduced using the Continuum and Line Anal- 
ysis Single dish Software (CLASS, a GILDAS software 1 ). In 
some of the sources, the CO mesospheric emission line was 
detected close to the cloud emission, which peaks at ui sr ^ 
+6km.s~', on the red side of the spectrum. For all our obser- 
vations we have adjusted a gaussian to the CO mesospheric line 
profile and subtracted it out. The CO mesopheric line profile 
is typically a few K bright, with a linewith of about 1 km.s -1 
(CHECK), much less than the velocity range of the cloud emis- 
sion and the outflow wing emission. When the outflow feature 
is on the blueshifted part of the line spectrum, the mesospheric 
CO is absolutely harmless. When the observed outflow wing is 
on the redshifted side, we checked that the CO mesospheric line 
is much narrower than the wing velocity range, hence does not 
hamper the detection of the latter. 
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Fig. 1. A Map showing the positions of the 1 1 spectra obtained for 
FU Tau. The CO(2-l) line is shown here. The scale of the grid is 0-15 
kms~' and -0.6-6 K in x and y respectively. For all spectra except the 
(0,-20) and (+20,+20) positions, the rms is less than 40 mK. For the 
(0,-20) and (+20,-1-20) positions the noise is ~ 80 mK. 



3. Results and Discussion 

3.1. Outflow signature 

In Fig. 2 a magnified view of the wings of the CO(2-l) emis- 
sion line at each point in the map of FU Tau (Fig. 1) is shown. 
The central (0,0) position spectrum is repeated as a dark solid 
line in all the plots, superimposed on the color lines observed at 
the other positions. As the beam is » 11" wide it encompasses 
both FU Tau A and FU Tau B. The bulk of the CO(2-l) emission 
comes from the cloud and is centered on the cloud velocity at 

1 http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. Montage of some of the CO(2-l) emission line profiles ob- 
served for FU Tau shown with a zoom on the region of the wings. The 
central (0,0) position spectrum is repeated as a dark solid line in all 
the plots, superimposed on the lines observed at the other positions, in 
color. We detect excess emission in the blue-wing between 3 kms -1 and 
5 kms" 1 , at the (0, 20), (10, 20) and (0, 30) positions. This points to a 
molecular outflow driven by the FU Tau system. 



6 km.s" 1 , and we search for outflow signature from variations in 
the CO(2-l) wing emission with respect to the emission on the 
central source. The outflow emission which is shifted in velocity 
with respect to the cloud is much fainter than the cloud emis- 
sion and therefore it will be detected in the wings of the CO(2-l) 
emission line. Fig. 2 shows that a blue component develops in 
the wing as we move towards the north, with an excess wing 
emission seen at the (0, 20), (10, 20) and (0, 30) positions be- 
tween a velocity of 3 kms -1 and 5 kms -1 . There is also a hint 
of blue-shifted excess emission at the (0, 10) position, and a red 
wing appears in the 8-10 km.s" 1 range in the (0,-20) position 
spectrum. The detection of this excess emission strongly points 
to a molecular outflow driven by FU Tau. In Fig. 3 we present 
an integrated intensity map of the blue-shifted wing emission in 
the velocity range 3 kms -1 to 6 kms -1 . The black squares mark 
the positions at which data was collected and we have super- 
imposed an optical image (WFCAM/UKIRT) of the FU Tau bi- 
nary taken from Luhman et al. (2009). The detection of the out- 
flow in the form of a blue-shifted lobe towards the north-east is 
clear. Without further data, we estimate a PA of ~ 20° for the 
outflow axis. 



3.2. Outflow Parameters 

Following Bachiller et al. (1990) we compute the CO column 
density in the blue lobe of the outflow with the following equa- 



PA i N(CO) M(H 2 ) M(H 2 ) 

O (°) (cm -2 ) (M ) (M Q /yr) 

20 60 3.6 10 16 4 ± 1 10~ 6 5 ± 1 10" 
Table 1. FU Tau outflow parameters 



tion. 

(CO) (cm -2) = 1.06 10 13 r mb exp(^) f T R (2 - l)dv (1) 

'mb J 

Adopting a gas excitation temperature T m \> ~ 15 K, with an 
H2/CO ratio of 10 4 and the results of Figure 3, we compute a 
mass in the blue lobe of the flow of M B (H 2 ) = 4 ± 0.8 10" 6 M Q . 

If we suppose that the momentum of the underlying jet has 
entirely been transferred to the molecular component that we ob- 
serve today, we can write : 

M B (H 2 ) < V max > = M jet < V jet > r dyn (2) 

We measure < V max >= 3 kms , and we take a canonical value 
< Vjet >= 100 kms -1 ; together with T dyn as 200 yr (see § 3.3, 
first paragraph), we obtain a mass loss rate for the blue-shifted 
lobe of M out = 6+1.3 10~ 10 M o /yr. This value can be modi- 
fied by various factors. For instance, we can adopt a correction 
factor to take into account the fact that the jet might have been 
episodic in the past. (Phan-Bao et al. 2011) use a factor of 10 
for this purpose. Also, the excitation temperature is uncertain 
although the T m ^ exp(16.5/r m b) factor does not vary much over 
7*21 = 10 - 25 K range. We could also take into account extinc- 
tion effects and the fact that we are only measuring half of the 
flow emission. Thus this value must be taken as a first order esti- 
mation of the outflow rate and most probably underestimates the 
rate of the underlying jet. The outflow parameters are summa- 
rized in Table 1 . 

3.3. Outflow powering source 

Although we cannot disentangle FU Tau A from FU Tau B it is 
most likely that the outflow is driven by the primary as forbidden 
emission associated with the primary has already been detected 
(Stelzer et al. 2010). Thus for the rest of the discussion we will 
assume that FU Tau A is the driver of the outflow. The peak of 
the blue lobe is measured at ~ 20" from the central source, pro- 
jected on the plane of the sky. Adopting a projection angle of 60° 
(Stelzer et al. 2013), the linear distance is « 50", corresponding 
to 7000 AU at the distance of Taurus (140 pc). At 100 kms -1 , this 
yields a dynamical age Td yn * 200 yrs for the observed outflow 
event. 

Previous to the results presented here, ISO-Oph 102 and 
MHO 5 with masses of 60 M Jup and 90 M Jup were the lowest 
mass objects for which molecular outflows were detected (Phan- 
Bao et al. 2008, 201 1). The outflow mass and mass outflow rate 
were estimated at M out = 1.6 x 10~ 4 M Q ; M out = 1.4 x 10~ 9 
M Q /yr and M out = 7.0 x 10~ 5 M ; M out = 9.0 x 10" 10 M /yr re- 
spectively. Thus our estimates of the mass and mass outflow rate 
of the FU Tau outflow agree with previous results and are in line 
with the fact that FU Tau A has the smallest mass of the three 
objects. The derived values of M out are also consistent with M out 
measured for the optical components of BD outflows (Whelan 
et al. 2009). For ISO-Oph 102 the outflow rate in the molecular 
component was found to be slightly higher than the the optical 
component. However, it is reasonable that M out for a molecular 
outflow could be greater than the outflow rate in the underly- 
ing jet. Assuming that the jet is powering the molecular outflow 
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(Downes & Cabrit 2007), the mass outflow rate of the molecu- 
lar component will grow with time as the jet transfers increasing 
amounts of energy and momentum. The size of the outflow of 
~ 20 " is compatible with the ISO-Oph 102 and MHO 5 flows 
and with observations of BD optical outflows. The agreement 
between the scale of the molecular and optical components is 
important if one is to accept that the jet drives the molecular 
flow. Finally we compare the mass outflow rate of the CO out- 
flow from FU Tau with the derived mass accretion rate and find 
the ratio of mass outflow to mass accretion, M out / M. dcc = 0.8 to 
1.7. M out I M acc measured for other BDs and VLMSs has also 
been found to be high compared to T Tauri stars where it is mea- 
sured at 1-10% (Bacciotti et al. 2011; Whelan et al. 2009). We 
stress that such a high value of M out / M acc cannot be used as a 
clue that this outflow removes or has removed a large fraction of 
the central object' mass (as in the Machida et al. (2009) model), 
because the observed FU Tau molecular outflow results from the 
entrainment by an underlying jet, and concerns an ejected mass 
several orders of magnitude smaller than in the (Machida et al. 
2009) model. On the same hand, another explanation for the for 
M out /M. d cc ratio « 1 in the current known series of BD sources 
could be due to an observational bias, because the first currently 
available observations are only sensitive to the most extreme jets 
in brown dwarfs. If this is true, the current ratios should prove 
much higher than the (as-yet unobserved) mean in BDs. More 
sensitive observations are thus needed to solve this issue. 




Fig. 3. An Integrated map of the blue-shifted wing emission in the 
velocity range 3-5 kms -1 . The LSR velocity of the BD is +6 kms" 1 . 
The black squares mark the positions at which observations we made. 
The level are 0.1, 0.15, 0.2, and 0.25 K.km.s -1 . Clearly we see a blue 
outflow lobe at a PA of ~ 20°. We have superimposed a UKIDSS K- 
band image taken from Luhman et al. (2009) of the BD binary at the 
same scale and assuming that FU Tau-A is at the central position. 



4. Summary 

The discovery of a molecular outflow driven by FU Tau A adds 
significantly to the interesting properties of this source and its 
binary companion. The FU Tau binary has a large separation 
compared to other binary systems and is thought to have formed 
in relative isolation. Both components harbor Class II accretion 
disks and FU Tau A is somewhat over-luminous for its spectral 
type. The fact that FU Tau A is driving an outflow demonstrates 
that despite having unusual characteristics it still exhibits prop- 
erties which are strongly linked to the formation of low mass 
protostars. The mass, scale and mass outflow rate we measure 
for the FU Tau A CO outflow agrees with previous observations 
of BD molecular outflows. While this result is a further impor- 
tant piece of evidence linking how BDs form to low mass star 
formation the derived ratio of mass outflow to accretion rates is 
much higher than what is observed in low mass protostars and 
in particular the T Tauri stars. For other BDs the two rates have 
been found to be comparable thus these new results for the FU 
Tau system support other studies of BD outflow activity (Bac- 
ciotti et al. 2011; Whelan et al. 2009). M out /M acc ratio « 1 
in the current known series of BD sources could be due to an 
observational bias, because the first currently available observa- 
tions could be only sensitive to the most extreme jets in brown 
dwarfs. More observations are needed to solve this issue. FU 
Tau is an excellent candidate for follow-up observations with 
sub-millimeter interferometer such as the Plateau de Bure in- 
terferometry or the Sub-miillimeter Array (SMA). With higher 
angular resolution observations we will be able to fully resolve 
the outflow, search for a red-shifted lobe and confirm whether 
FU Tau A is the driving source of the flow. 
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